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Abstract. We have used the Hubble Space Telescope to measure the motion in the sky and compute the galacto- 
centric orbit of the black hole X-ray binary GRO J1655-40. The system moves with a runaway space velocity of 
112 ± 18 km s“^ in a highly eccentric (e = 0.34 ± 0.05) orbit. The black hole was formed in the disk at a distance 
greater than 3 kpc from the Galactic centre and must have been shot to such eccentric orbit by the explosion of 
the progenitor star. The runaway linear momentum and kinetic energy of this black hole binary are comparable 
to those of solitary neutron stars and millisecond pulsars. GRO J1655-40 is the first black hole for which there is 
evidence for a runaway motion imparted by a natal kick in a supernova explosion. 
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1. Introduction 


Neutron stars are known to have large transverse motions 
on the plane of the sky which are believed to result from 
natal kicks imparted by supernova explosions. Energetic 
explosions have also been invoked in models of the core 
collapse of massive stars onto black holes. However, there 
have been few observations that constrain models of the 
physical processes by which stellar-mass black holes are 
formed. The measurement of a large radial velocity for the 
centre of mass of the black hole X-ray binary GRO JI655- 


40 (Orosz & Bailyn, 1997; Shahbaz et ah, 1999), together 
with chemical elements found by Israelian et al. (1999) 
in the surface of the donor star, provided observational 
support to the idea that black holes -as neutron stars- 
may form in supernova explosions that impart strong natal 
kicks to the collapsed objects. 

If a black hole is accompanied by a mass-donor star, it 
is possible to determine the radial velocity, proper motion, 
and distance of the system, from which one can derive the 
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space velocity, track the path to the site of birth, and con¬ 
strain the strength of the natal kick. Presently, the most 
accurate proper motions of X-ray binaries are obtained 
following at radio wavelengths with Very Long Baseline 
Interferometry (VLBI) the motion in the sky of the asso¬ 
ciated compact microquasar jets, as done recently for the 


halo black hole binary XTE J1118-1-480 ( Mirabel et al. 
[2001 ). This has not been possible for GRO J1655-40 be¬ 


cause there is no VLBI calibrator nearby, and in the last 
years the radio counterpart faded away below detection 
limit. In this context, we carried out relative astrometry 
of the secondary star using optical images obtained with 
the Hubble Space Telescope 6.3 years apart. Here we re¬ 
port the proper motion of GRO J1655-40 which together 
with the radial velocity allows us -taking into account the 
uncertainties in the distance- to determine the parameters 
of its runaway kinematics and galactocentric orbit. 


2. The proper motion of GRO J1655-40 

On April 1996 the field of GRO J1655-40 was observed 
for the first time with the WFPG2 on the Hubble Space 


Telescope (Tavani et al., 1996). In quiescence, the sec¬ 
ondary star has an apparent magnitude of my ~ 17.2. 
The target had been placed at the center of the Planetary 
Camera (PC) chip, with a pixel size of (X'OdS. Observations 

































2 


I.F. Mirabel et al.: The runaway black hole GRO J1655-40 



Fig. 1. Position of GRO J1655-40 on a R band image from the Digitized Palomar Observatory Sky 
Survey II (POSS II), in Galactic coordinates. The arrow shows the direction of the motion at a rate of 5.2 ± 
0.5 mas yr“^ measured with the Hubble Space Telescope. Most of the stars around I = 345.44°, b = + 2.43° belong 
to the open cluster NGC 6242 which is at a distance of 1 ± 0.1 kpc from the sun. Because of the uncertainty in the 
distance to GRO J1655-40 the association with the cluster cannot be assessed. 


were obtained through the 675W filter (A = 6717 A; 
AA = 736 A). A total of 16 exposures of 40 s each were ac¬ 
quired on April 26*^ and the same sequence was repeated 
on April 27*^ but with a different telescope roll angle. Both 
sets of exposures were taken in groups of 4 and each group 
was dithered by a few pixels in both Right Ascension and 
Declination with respect to the others. 


On June 20*^ 2001 a total of 18 exposures of 40 s each 
were acquired by us, with the same observational set-up 
as in April 1996, but without any dithering between single 
exposures. After a pre-reduction with the standard HST 
pipeline reduction (debiasing, dark removal, flatfielding), 
groups of well aligned exposures (i.e. with relative shift 
smaller than 0.01 pixel) were combined with a median 
filter, stacked and cosmic ray hits filtered out. 


We finally ended up with 4 images for each of the 
two April 1996 observations and one image for the June 
2001 one. All the 1996 images were then registered on 
the 2001 one, previously aligned in Right Ascension and 
Declination, by htting the coordinate transformation be¬ 
tween grids of reference objects. The correction for the 
WFPG2 geometrical distortions, optimized for the used 
filter (Trauger et al., 1995) was taken into account. The 
whole procedure was iterated until the transformation 
residuals were all below ~ 1.5 tr. We finally computed the 
relative displacements of our target for each of the 8, inde¬ 
pendent, couples of 1996/2001 images and we averaged the 
results. The computed proper motion is /io, = —3.3 ± 0.5 


..-1 


mas yr^" and fis = —4.0 ± 0.4 mas yr“^, corresponding 
to an overall yearly displacement /r = 5.2 ± 0.5 mas yr“^ 
along a Position Angle of 220° ±5°. At a distance D(kpc) 
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this proper motion corresponds to a transverse velocity on 
the plane of the sky of (25±3) D(kpc) km s“^. In Figure |] 
we show the path of the black hole binary on a -R-band 
image of the Digitized Palomar Observatory Sky Survey 
II (POSS II). NGC 6242 is a well studied open cluster at a 
distance of 1020 ± 100 pc from the Sun (plushkova et al 


1997). Figure shows that GRO J1655-40 is close to the 


two equations (Mirabel & Rodriguez, 1994): 


A^r,a — 


l3sin{0) 


1 ± l3cos(9) D ’ 


( 1 ) 


where and fj,a are the proper motions of the receding 
and approaching jets. In this system of two equations there 
are three unknowns: the angle with the line of sight of 
the jet axis 9, the velocity of the jets f3 = ^, and the 
distance D. We point out that using solely the observa tions 
at radio wavelengths by Hjellming fc Rupen ( 1995 ) it is 
not possible to solve these equations, unless on e assumes 
a value fo r one o f the three variables 9, or D. Hjellming 
fc Rupee ( 1995 ) assumed 9 = 85° ± 2°, from which one 
de rives /3 = ^ = 0.92 and D = 3.2 kpc. As already noticed 
by Orosz fc Bailyn ( 1997 ), in this case the axis of the jet 
and the axis of the orbital plane differ by ~ 15° ± 2°. 

We point out that the assumption that the jet axis is 
parallel to the axis of the orbital plane is equally consistent 
with the observations at radio wavelengths. The jet axis 
and the orbital plane must be coupled, since the period of 
rotation of the jets about the jet axis (Hjellming fc Rupen 


boundary of a dark cloud (cataloged as DC344.9+2.6 in 
Hartley et al.^ |1986 ), which could explain the relative large 
reddening of the secondary star in GRO J1655-40. 

3. The distance of GRO J1655-40 

In order to use our proper motion measurement to con¬ 
strain the space velocity of the black hole, we need to 
know its distance. Unfortunately, as for most X-ray bina¬ 
ries, the distance to GRO J1655-40 is rather uncertain. 
In the following we discuss the observational constrains to 
the distance of this source. 

3.1. The relativistic distance 

It is widely believed that the distance of GRO J1655-40 
was well determined solely from the kinematics of the two 
sided radio jets. This is not true. The relativistic time 
delay of the motion of the ejecta in the sky is given by the 


In summary, from the data at radio wavelengths and 
the system of two equations (|^) with three unknowns 
one can only derive with certainty a relativistic up¬ 
per limit ( Mirabel fc Rodrigue4 |1999| ) given by D < 
c/(/rrA*a)“^''^ = 3.5 kpc. 


3.2. The distance and the interstellar matter along the 
line of sight 

A distance for GRO J1655-40 was propose d on the basis 
of optical (Bailyn et al.|, |1995 ) and X-ray ( [Greiner et al. 


1995; Ueda et al., 1998a|,b|) measurements of the column of 


interstellar matter in the line of sight, under the assump¬ 
tion that the absorbing material is distributed homoge¬ 
neously between the source and the observer. However, 
GRO J1655-40 is at relatively high Galactic latitude 
(Galactic longitude and latitude I = 345.0°, b = -|-2.2°) 
in the Scorpius region of the sky which contains rather 
clumpy optical dark clouds in the foreground (see Figure 
that have 60^m and lOO^m IRAS counterparts of dust 
emission. From a study of the reddening undergone by the 
stars that are at distances between 700 pc and 1900 pc, it 
is known that most of the reddening in this region of the 
sky occurs in the local arm within 700 pc from the Sun 
(Crawford et al., 19891 ). 


On the other hand, a kinematic distance was proposed 
from the radial velocity of absorption features in the HI 
A21cm line spectrum (Tingay et ah, 19951 ). However, it is 
known that in this region of the sky at distances < 1900 
pc there are clouds with anomalous velocities of up to 


—50 km s ^ (Crawford et al, 1989| ). Therefore, it is not 


possible to derive the distance of GRO J1655-40 only from 
the column density and/or kinematics of the interstellar 
matter in the line of sight. 

3.3. Constrains on the distance from the properties of 
the secondary star 

The main argument in favor of the canonical distance of 
^3.2 kpc is based on the flux, color and size of the sec¬ 
ondary. The radius is inferred from the photometric light 
curves which provide evidence that the secondary fills its 
Roche lobe (prosz fc Bailyn , 1997; shahbaz et al, 1999 


Soria et al, ^000| ; van der Hooft et al, 1997, 1998; Phillip; 


3t al, 19991). From the optical spectrum (Orosz & Bailyn 


1997) and interstellar absorption (Horne et al, 1996]) a 


1995) is -within the uncertainties-, the same as the 2.6 day 


orbital period of the binary ( prosz fc Bailyn , 1997). If one 
assumes that the twin jets are perpendicular to the orbital 
plane of the binary, from pr = 45 mas/day and Pa = 54 
mas/day ([Hjellming fc Rupen , 1995), 9 = 70.2° ± 1.9° 
( preene et al , 2001) result (3 = ^ = 0.27 ± 0.03 (where 
V is the velocity of the jets and c the speed of light) and 
a distance D = 893 ± 100 pc. Under this assumption, the 
distance would be a factor 3.5 closer than commonly as¬ 
sumed, the jet velocity would be similar to that in SS433, 
and GRO J1655-40 would not be a superluminal source. 


temperature is derived, which together with the radius 
pr ovides an intrinsic lum i nosity . In the most recent model 
by Beer fc Podsiadlowski ( 2002[ ) the secondary would have 
a luminosity of 21 ± 6.0 Lq, which is consistent with a 
distance >2 kpc. These authors estimate masses for the 
black hole VIbh = 5.4 ± 0.3 Mq and for the donor star 
M* = 1.45 ±0.35 Mq. 

We point out that the secondary star in quiescence has 
an apparent magnitude my = 17.12, it has been classified 


as an F3 IV-F6 IV sub-giant (Orosz & Bailyn, 1997), and 
along the line of sight there is an interstellar absorption 
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Ay = 3.1 X E(B — V) = 4.03 mag ( Horne et al. , 1996| ). A 
sub-giant star of this spectral type has a mean intrinsic 
magnitude My ~ 3.2±0.2 mag (Popper, 198C; Aller et al 


1982), and it would be at a distance D = 950 ± 150 pc. 


Alternatively, if the secondary were a main sequence star 


of spectral type F5 V star (Regos et al,, 1998), for an 


absorption Ay in the range of 3.3-4.3 mag the distance 
would be in the range of 800-1250 pc. However, it may be 


D [kpc] 

0.9 

3.2 

Vrun [km 8“-^] 

130 

93 

e 

0.39 

0.29 

^max [kpc] 

0.05 

0.15 

Tmax [kpc] 

13.8 

7.2 

Tmin [kpc] 

6.0 

3.9 

p [Mq km s“^] 

430 

637 

Tkin [erg] 

5.6xl0'‘" 

5.9x10'*^ 


incorre 
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Fol 
Podsia 
luminop: 
kpc wi 
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star ( 
secondk: 
mass 0 


ct to attribute the absolute magnitudes of isolated 
) secondary stars in X-ray binaries with the same 
1 type. 

owing kind communications by Beer & 
llowski (private communication), the lower 
,ity and temperature implied by a distance of ^1 
require a higher mass ratio with much reduced 
for the compact object (^3.2 Mq) and secondary 
'0.1 Mq). But the absorption spectrum of the 
ry seems to rule out an M type star of such low 
r a K-tvpe star with mass 0.6 <Mf 7 ^. The spectra 


culd 


plane, and 


Table 1. Parameters of the runaway motion and 
galactocentric orbit of GRO J1655-40, for two ex¬ 
treme val ues of the distance. Vrun is the runaway ve¬ 
locity in tliree dimensions after correction for differential 
Galactic rotation, e is the eccentricity of the galactocen¬ 
tric orbit, ^max is the maximal height above the Galactic 


and 


are the maximal and minimal 


galactocentric distances. The runaway linear momentum 
p and kinetic energy Tkin of the binary system are com¬ 
puted according to the masses given in the text. 


of stars with Te//<5000 K have different signatures, such 
as strong molecular bands. This was the case in GRS 


1915±105 where K band spectroscopy (Greiner et al 


2001 ) revealed CO molecular bands rendering invalid the 
cla.ssifipation of the donor 3,s a, mfl.iri sequence star 


Soderblom (1987)’s equations of transformation, and as¬ 
suming the sun moves (Uo,V 0 ,Wo) = (9, 12, 7) km s~^ 
relative to the local standard of rest (l.s.r.) (Mihalas fc 


We point out that there has been arialngniis uncer- Binney , 1981 ). Two different orbits were computed for two 


tainties about the nature of the secondary in the X-ray 
binary LMX-3; the proposition that it is a main sequence 


star ( 

Cowley et al, 

00 

an der Klis et al. 

1985|) has 

recently been challenged by 

Soria et al. (^001 

) who ar- 


gue that it is a sub-giant. Furthermore, it is not known 
what could be in GRO J1655-40 the contribution to the 
optical flux from: 1) the accretion disk detected in the X- 
rays at times when it was believed that the source was in 
optical quiescence ( [Garcia et ah , 200lD , and 2) possible 
non-thermal processes (e.g. synchrotron jets) that may be 
associated with the polarization observed in the optical 
flux (Gliozzi et al, 1998| ). In this context, for the scope of 
the present study we leave as an open question the issue 
on the distance of GRO J1655-40. _ 


extreme values of the distance: for D = 0.9 kpc we obtain 
(U,V,W) = (-133 ± 2, 27 ± 2,1 ± 3) and for D = 3.2 kpc, 
(U,V,W) = (-121 ± 18, -33 ± 8 ,3 ± 8 ). These two sets of 
values are rather different from the mean values that char¬ 
acterize the kinematics of stars that belo ng to the halo , 
and t he thin and thick disk of the Galaxy ( jChiba &: Beers , 
^OOPj) . The runaway velocities Vrun were computed for the 


two possible extreme distances of 0.9 kpc and 3.2 kpc (see 
Table 1), after subtracting the Galactic differential rota¬ 
tion given by the model for the corresponding position of 
the source in the Galactic disk. 

The runaway linear momentum p and kinetic energy 
Tfein were computed assuming Mbh = 5.4 ± 0.3 Mq and 
M, = 1.45 ± 0.3 Mq for a distance D = 3.2 kpc (Beer & 
Podsiadlowski, 2002| ), and M^h = 3.2 Mq and M* = 0.1 


4. Space velocity and galactocentric orbit 

As shown below, the parameters of the runaway motion 
and galactocentric orbit of GRO J1655-40 are essentially 
the same for distances in the range of 0.9-3.2 kpc. The 
tangential velocity is 


Vt = (25 ± 3) Dkpckms 


( 2 ) 


where Dkpc is the distance in kpc. The radial velocity with 
respect to the Sun is -142.4±1.5 km s“^ (prosz &: Bailyn, 


1997; Shahbaz et al, 1999| ). 


Using values of the position, distance, proper motion 
and radial velocity, the Galactic orbit of GRO J1655-40 
can be computed using a Galactic gravitational potential 
model (Johnston et al, 1995] ). The velocity components U, 
V, and W directed to the Galactic centre, rotation direc¬ 


tion, and north Galactic pole are derived using Johnson & 


Mq for a distance D = 0.9 kpc (Beer & Podsiadlowski, 
private communication). 

The parameters of the runaway motion and galacto¬ 
centric orbit of GRO J1655-40 are given in Table 1. Jn 
Figure || are represented the Galactocentric orbits. For 
a given Galactic potential the orbital parameters do not 
change as a function of time. The selection of different 
Galactic potentials from the current models introduce a 
scatter in the values of the orbital parameters smaller than 
10%. Jn fact, the errors in the parameters listed in Table 1 
are largely dominated by the uncertainty in the distance. 
For the range of distances 0.9-3.2 kpc, the galactocen¬ 
tric orbit is highly eccentric (e = 0.29-0.39), the source 
always moves within a maximum height of 150 pc from 
the Galactic plane, and since the minimum perigalactic 
distance is 3.9 kpc it never reaches the Galactic bulge. 

Clearly, GRO J1655-40 does not move in a halo orbit 
as XTE J1118±480 ( [Mirabel et ^ p00T[ ). GRO J1655-40 
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Fig. 2. Galactocentric orbits of GRO J1655-40 viewed from above the Galactic plane. The orbital plane 
is almost parallel to the Galactic disk and the source never reaches heights greater than 150 pc. On the left is shown 
the orbit obtained for an heliocentric distance D = 0.9 kpc, and on the right, for D = 3.2 kpc. GRO J1655-40 never 
penetrates the Galactic bulge. Plotted in dashed blue line is 1 Gyr of backwards integration, and in red the past 230 
Myrs. 


must have been born in the Galactic plane at a galacto¬ 
centric distance > 3.8 kpc. The runaway linear momentum 
and kinetic energy of the binary system are comparable 
t o those of solitary n eutron stars and millisecond pulsars 
( fPoscano et ah , 1999 ). 


5. Conclusion 

From the proper motion and radial velocity of GRO J1655- 

40 -irrespective of the uncertainties in the distance- we 

conclude the following: 

1. The X-ray binary has a runaway velocity of 112±18 km 
s“^, probably imparted by a natal explosion during the 
formation of the black hole. 

2. The linear momentum of the binary system is 538 ± 
100 Mq km s“^, which is comparable to that found 
in solitary neutron stars and millisecond pulsars. This 
suggests that the relatively low-mass black hole (< 6 
Mq) in GRO J1655-40 may have been formed through 
a neutron star phase. 

3. The kinetic energy of the system is (5.8±0.4) x 10^^ 
ergs, namely, ~6 x 10“"* that of a typical supernova. 

4. The galactocentric orbital plane is almost parallel to 
the Galactic disk and the Galactic pericentre > 3 kpc, 
which indicates that the black hole in GRO J1655-40 
is the remnant of a massive star born in the Galactic 
disk. 

5. Contrary to the binary XTE J1118-I-480 which harbors 
a black hole born in the halo ( [Mirabel et al. , 2001), in 
the case of GRO J1655-40 there is unambigous kine- 
matical evidence that the black hole was born in the 


disk and received a natal kick, most like ly from a su¬ 
pernova explosion (Israelian et al., 1999|). 


Acknowledgements. We thank J. Casares, P. Podsiadlowski, 
M.E. Beer, J.-P. Lasota, S. Corbel, P. Benaglia, A. Piatti, 
and E. Gourgoulhon for help and discussions on different as¬ 
pects of this work. J.A.C. and I.F.M. are members of the 
Consejo Nacional de Investigaciones Cientificas y Tecnicas 
of Argentina, and I.R. is a fellow of the Conselho Nacional 
de Desenvolvimento Ciem'fico e Tecnologico of Brazil. I.F.M. 
acknowledges support from PIP 0049/98 and Fundacion 
Antorchas. 


References 

After, L. H., Appenzeller, L, Baschek, B., et ah, eds. 1982, 
Landolt-Bornstein: Numerical Data and Functional 
Relationships in Science and Technology - New Series 
” Group 6 Astronomy and Astrophysics ” Vol.2 
Schaifers/Voigt: Astronomy and Astrophysics” Stars 
and Star Clusters 

Bailyn, C. D., Orosz, J. A., Girard, T. M., et al. 1995, 
Nature, 374, 701 

Beer, M. E. & Podsiadlowski, P. 2002, MNRAS, 331, 351 

Chiba, M. & Beers, T. C. 2000, AJ, 119, 2843 

Cowley, A. P., Crampton, D., Hutchings, J. B., Remillard, 
R., & Penfold, J. E. 1983, ApJ, 272, 118 

Crawford, I. A., Barlow, M. J., & Blades, J. C. 1989, ApJ, 
336, 212 

Garcia, M. R., McClintock, J. E., Narayan, R., et al. 2001, 
ApJ, 553, L47 



























6 


I.F. Mirabel et al.: The runaway black hole GRO J1655-40 


Gliozzi, M., Bodo, G., Ghisellini, G., Scaltriti, F., & 
Trussoni, E. 1998, A&A, 337, L39 
Glushkova, E. V., Zabolotskikh, M. V., Rastorguev, A. S., 
Uglova, 1. M., & Eedorova, A. A. 1997, Astronomy 
Letters, 23, 71 

Greene, J., Bailyn, G. D., & Grosz, J. A. 2001, ApJ, 554, 
1290 

Greiner, J., Guby, J. G., McGaughrean, M. J., Gastro- 
Tirado, A. J., & Mennickent, R. E. 2001, A&A, 373, 
L37 

Greiner, J., Predehl, R, & Pohl, M. 1995, A&A, 297, L67 
Hartley, M., Tritton, S. B., Manchester, R. N., Smith, 
R. M., & Goss, W. M. 1986, A&AS, 63, 27 
Hjellming, R. M. & Rupen, M. P. 1995, Nature, 375, 464 
Horne, K., Harlaftis, E. T., Baptista, R., et al. 1996, 
lAU Circ., 6406, 2 

Israelian, G., Rebolo, R., Basri, G., Gasares, J., & Martin, 
E. L. 1999, Nature, 401, 142 

Johnson, D. R. H. & Soderblom, D. R. 1987, AJ, 93, 864 
Johnston, K. V., Spergel, D. N., & Hernquist, L. 1995, 
ApJ, 451, 598 

Mihalas, D. & Binney, J. 1981, Galactic astronomy: 
Structure and kinematics /2nd edition/ (San Francisco, 
GA, W. H. Freeman and Go., 1981. 608 p.) 

Mirabel, I. F., Dhawan, V., Mignani, R. P., Rodrigues, I., 
& Guglielmetti, F. 2001, Nature, 413, 139 
Mirabel, I. F. & Rodriguez, L. F. 1994, Nature, 371, 46 
—. 1999, ARA&A, 37, 409 
Grosz, J. A. & Bailyn, G. D. 1997, ApJ, 477, 876 
Phillips, S. N., Shahbaz, T., & Podsiadlowski, P. 1999, 
MNRAS, 304, 839 

Popper, D. M. 1980, ARA&A, 18, 115 
Regos, E. ., Tout, G. A., & Wickramasinghe, D. 1998, ApJ, 
509, 362 

Shahbaz, T., van der Hooft, F., Gasares, J., Gharles, P. A., 
& van Paradijs, J. 1999, MNRAS, 306, 89 
Soria, R., Wu, K., & Hunstead, R. W. 2000, ApJ, 539, 445 
Soria, R., Wu, K., Page, M. J., & Sakelliou, I. 2001, A&A, 
365, L273 

Tavani, M., Fruchter, A., Zhang, S. N., et al. 1996, ApJ, 
473, L103 

Tingay, S. J., Jauncey, D. L., Preston, R. A., et al. 1995, 
Nature, 374, 141 

Toscano, M., Sandhu, J. S., Bailes, M., et al. 1999, 
MNRAS, 307, 925 

Trauger, J. T., Vaughan, A. H., Evans, R. W., & 
Moody, D. G. 1995, in Galibrating Hubble Space 
Telescope. Post Servicing Mission. Proceedings of a 
Workshop held at the Space Telescope Science Institute, 
in Baltimore, Maryland, May 15-17, 1995. Editor(s), 
Anuradha Koratkar, Glaus Leitherer; Publisher, Space 
Telescope Science Institute, Baltimore, Maryland, 1995. 
LG #: QB500.268 .G34 1995. ISBN #: NGNE., p.379, 
379-386 

Ueda, Y., Inoue, H., Tanaka, Y., et al. 1998a, ApJ, 492, 
782 

—. 1998b, ApJ, 500, 1069 


van der Hooft, E., Groot, P. J., Shahbaz, T., et al. 1997, 
MNRAS, 286, L43 

van der Hooft, F., Heemskerk, M. H. M., Alberts, F., & 
van Paradijs, J. 1998, A&A, 329, 538 
van der Klis, M., Glausen, J. V., Jensen, K., Tjemkes, S., 
& van Paradijs, J. 1985, A&A, 151, 322 



